Superconductor-to-normal transition in finite nanowires 
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In this paper we discuss the interplay of quantum fluctuations and dissipation in uniform super- 
conducting nanowires. We consider a phenomenological model with superconducting and normal 
components, and a finite equilibration rate between these two-fluids. We find that phase-slip dipoles 
proliferate in the wire, and decouple the two-fluids within its bulk. This implies that the the normal 
fluid only couples to the superconductor fluid through the leads at the edges of the wire, and the 
local dissipation is unimportant. Therefore, while long wires have a superconductor-metal transition 
tuned by local properties of the superconducting fluid, short wires have a transition when the total 
resistance is iitotal = Rq = h/Ae 2 . 
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Quantum phase transitions have long been at the fore- 
front of condensed matter theory. Especially interesting 
are systems of reduced dimensionality and size, where 
fluctuations are enhanced, and ordering is illusive, and 
far from being expalined by mean field theory. Such sys- 
tems exhibit a surprising degree of universality; for in- 
stance, as observed in Refs. [1, 2], a mesoscopic Joseph- 
son junction shunted by a resistor R undergoes a (so- 
called Schmid) transition between a Coulomb-blockade 
(normal) and superconducting phase when the shunt re- 
sistor is R = Rq = h/Ae 2 = 6.A5kQ [3-5]. Fluctuations 
of the superconducting phase angle, i.e., phase-slips, in- 
duce this transition; they also control the onset of su- 
perconductivity in long thin wires and Josephson junc- 
tion chains [6-10], where the competition between local 
charging energy, which creates phase-slips, and the su- 
perconducting stiffness tunes the transition [11-13]. 

Our focus is experiments on Mo^Gei\ (amorphous) 
nanowires as narrow as 5nm-15nm. Resistance vs. tem- 
perature curves showed a transition between supercon- 
ducting (resistance decreasing upon cooling) and normal, 
or weakly insulating, (resistance non-decreasing upon 
cooling) behavior. A first set of measurements on wires 
of various diameters, and lengths lOOnm < L < 200tito, 
showed a remarkable result: a transition when the to- 
tal resistance of the wire was Rq = h/Ae 2 = 6. A5kH, 
[14], as if the entire wire was a single shunted junc- 
tion. But the coherence length of MoGe is £ < IOtito, 
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FIG. 1: A two- fluid model of a superconducting nanowirc 
with a normal part, depicted as a separate region. In order 
for electrons to change from normal to superconducting, they 
need to pass through the conversion layer, which has a finite 
conductivity. Proliferated phase-slip dipoles inhibit two-fluid 
relaxation, and make the conversion layer insulating, which 
renders the normal part as an effective single shunt resistor. 



thus the wire should differ dramatically from a single 
junction. Indeed, later experiments on longer wires, 
200nm < L < lOOOnm, showed a weak transition that 
depended on the resistance per length or cross-section of 
the wires, i.e., on a local quantity, rather than the to- 
tal resistance [15]. Later experiments [16], could neither 
prove nor disprove the global nature of the transition in 
the shorter wires. 

In this paper, we describe nanowires using a two-fluid 
model, which assumes that Cooper pairs couple to a nor- 
mal electron fluid, which provides local dissipation (Fig. 
1). Remarkably, we find that at sufficiently low tem- 
perature, the normal and superconducting fluids within 
a continuous nanowire decouple due to quantum phase 
fluctuations, thus rendering the local dissipation unim- 
portant. Therefore, the superconducting degrees of free- 
dom can only couple to the dissipative normal fluid at 
the leads, on the edges of the wire, where they couple 
to its total normal-state resistance, i?totai- As a result, 
we show that indeed short wires may undergo a global 
dissipative Schmid transition tuned by the total wire re- 
sistance, when i?totai = Rq- By short, we mean wires 
with length L > £, but shorter than both the thermal 
length, L < Tic/T (with c the Mooij-Schon velocity, and 
T being the lowest temperature in the experiment) [17], 
and the 'quantum length' AqRqI p, (with p the specific 
resistance, and Ac the largest cross-section area where 
quantum phase slips are not competely suppressed) [18]. 
After our work was completed, this result was verified in 
MoGe wires with 50nrn < L < 300nm [19]. Below we 
derive the two-fluid model, show how phase-slip dipoles 
decouple the normal and superconducting fluids, and ap- 
ply to model to the case of a finite nanowire. 

The hint of a Schmid transition, the long resistive tails 
seen in experiments, and the strong disorder of the M oGe 
nanowires suggest the presence of local dissipation, which 
motivates the two-fluid model approach. We assume that 
charge can flow in the nanowires in two ways: as diffusive 
normal electrons with resistivity p - normal fluid - and 
as bosonic Cooper pairs - superfluid. The normal fluid 
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stems from strong disorder and phase- fluctuations, which 
suppress the proximity effect and possibly give rise to 
normal regions and a finite density of states for single 
electrons at the Fermi level. The two fluids can have 
a different chemical potential, and can exchange charge 
with a finite, bare, relaxation time, r r = T _1 , in a bulk 
system (see Fig. 1, and Fig. 2b for a discrete model). 
This is related to the branch imbalance relaxation time 
[20, 21], first measured by Clarke [22] in Sn wires. 

Before plunging to the analysis, note that earlier 
works on similiar models considered only the perfect 
normal-super fluid coupling, T = oo case, and found a 
superconducing-metal transition tuned by the resistance 
per length [11-13, 23], as did Refs. [24, 25]. Alternative 
approaches assumed external dissipation coupled to the 
leads but not to the bulk of the wire [26], or discussed 
the onset of superconducting correlations and neglected 
phase fluctuations [27]. 

Indeed, in our model as well, sufficiently long but fi- 
nite wires should exhibit a SC-normal crossover tuned by 
their cross-section area, which sets the bare fugacity of 
quantum phase slips [18], as well as their stiffness. But 
quantum fluctuations in the form of phase-slip dipoles, 
make the Cooper-pair to normal-electron conversion rate 
vanish at T = in the bulk of the nanowire: T — > 0. As 
claimed above, this leads to a true Schmid transition for 
short wires, which effecively become a short dissipation- 
less superconducting wire, shunted through the leads by 
the total normal-state resistance Rt- 

The crucial two-fluid decoupling is already evident in 
a simple two-junction system (Fig. 2a) [28, 29]. When 
r = (i.e. vanishing conversion resistance), the two junc- 
tions in the system are independent in the d.c. limit. 
Phase slips - events where the phase across a Josephson 
junction tunnels by 2ir - create a sudden voltage drop 
that opposes any supercurrent flowing, and thus induce 
dissipation. A Schmid transition occurs in each junction 
when Ri = Rq (i = 1,2). When r > 0, the two junc- 
tions become coupled, and phase-slips may form bound 
dipoles: simultaneous phase-slip and anti-phase-slip in 
the two junctions. Remarkably, dipoles do not destroy 
the coherence between the two leads, since they produce 
equal and opposite voltage drops. Nevertheless, as sin- 
gle phase slips block supercurrents across their Joseph- 
son junctions when they proliferate, dipoles block the 
normal-superfluid conversion channel: a conversion cur- 
rent 2i (Fig. 2a) flowing across r, with no lead-to-lead 
current, implies a current i on both junctions, but in op- 
posite directions, i couples directly to the voltage drop 
of the dipoles; when proliferated, they block this current 
mode, and thus decouple the normal and super fluids. 
Phase slip dipoles proliferate roughly when r > Rq. In 
this case a global Schmid transition takes place when 
R\+ R% = Rq. 

Next, we generalize the normal-super fluids decoupling 
to wires, first using a discrete model (Fig. lb), and 
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FIG. 2: (a) Single two-fluid grain (pink ellipse) between two 
leads: dipoles produce two opposing voltage spikes on the 
two junctions, which oppose super-to-normal conversion cur- 
rents leaving the superconducting part of the grain (bottom 
circle in ellipse) and entering the normal part (top circle). 
When dipoles proliferate, the superfluid-normal conversion 
resistance (and time) effectively diverges, r — » oo, and the 
normal and superfluid are completely decoupled, (b) We be- 
gin our study with a chain of two fluid grains, which is a 
discretized version of the nanowire in Fig. 1. Note that 
1/r = Ta x and R = a x p. 



then taking its continuum limit. Starting with an infi- 
nite chain of mesoscopic two- fluids grains (Fig. lb) [30], 
the low-energy action for the chain is given in terms of a 
2d gas of phase-slips, with interaction: 



P1P2 K log 



\Jx 2 jc? + ' 



+ ae -MM<j , £l 



(1) 



c = a x \JEjEc /h is the Mooij-Schon velocity [31], and 
a T = a x /c. pi = ± is the phase slip polarity. The first 
term is the usual isotropic interaction of a 1+1 XY model 
due to the plasmons in the Josephson junction array; 
K = 2-Ky/Ej/Ec, E c = (2e) 2 /C. The second term is 
due to the dissipative interaction: a = max{ ^= , }, 

/Tp, a x } is a new length- 



and \q = max{a x y/r/R = I/a 
scale that arises from the two-fluid finite relaxation time. 
As in the two-junction case, dipoles must be explicitly 
included in the low-energy description of this model [12, 
13, 30]. We denote the fugacity of single phase slips as 
C, and the fugacity of a dipole with moment n as r] n . For 
completeness, we quote here the explicit field theory for 
the infinite chain: 
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with c = a x \/EjEc /h~, and 8, ip mediating the plasmon 
and dissipative interactions, respectively. At high ener- 
gies r\ n ~ C 2 , and A n /, = f i+n - /,-. This is a rep- 
resentation dual to the SC phase representation, hence 
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FIG. 3: To describe a continuous wire, we take the limit a x — > 
while keeping the coherence length and size of phase slip, 
£, fixed. This implies that phase slips are spread over £/a x 
junctions, and can form dipoles with separation x < £. The 
voltage signs symbolize the voltage drop caused by a phase 
slip. 



exp(i^j + iOi) is the operator that creates a phase slip on 
junction i. 

It is useful to compare the relatively complicated inter- 
action between phase slips in an infinite chain, with that 
of phase slips in a single Josephson junction. In a single 
junction the interaction is: p\P2-^§- log fe. The Schmid 
transition, which marks phase-slip proliferation, occurs 
when the gain in entropy due to separating a phase-slip 
from an anti phase slip, S = log(a T T) equals the re- 
quired interaction energy, -^-log(a T T). Employing the 
same argument for the two-fluid Josephson chain yields 
the approximate SC-normal phase boundary: K + ~ 
4. This transition is essentially the 1+1 Kosterlitz- 
Thouless-Berezinski (KTB) transition of the Josephson 
junction array in accordance with Rcf. [11] (for a more re- 
fined analysis see Ref. [30]). But this argument, as well as 
the interactions in (1), ignores phase-slip dipoles. When 
dipoles are proliferated, the normal- and superfluids de- 
couple, and the superconducting part of the wire exhibits 
a SC-normal KTB transition when K ~ 4. Phase-slip 
dipoles, we find, proliferate when: 
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The left-hand side is the strength of dipole interaction, 
which consists of the self interaction of the slip and anti 
slip and also their mutual interaction. a x is the distance 
between grains in the model. 

In the continuum limit, dipoles always proliferate and 
cut off the superfluid-normal conversion. The contin- 
uum limit implies a x — -> 0; but this makes a Joseph- 
son junction (and therefore also a phase slip on a junc- 
tion) shrink to length zero. But a phase slip occur- 
ring on physical nanowires has a characteristic length 
£ (coherence length). To reconcile this we allow phase 
slips to smoothly spread over ~ i/a x junctions [30]. 
Technically, we transform the zeta term in Eq. (2) as 
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where 



f(r) is a smooth, normalized, function centered around 0, 
with width £. We similarly treat the dipole r\ n terms. The 
smearing reflects that in nanowires, phase slips can have 
an almost arbitrary overlap, A < £ (Fig. 3), with other 
phase slips. The continuum generalization of Eq. (3) 



is that dipoles proliferate when: ^ ^ max {^ Q ^ 

(where i? 5 = R to tn\£,/L = p£). Thus for any r 
there is a separation A c below which dipoles prolifer- 
ate. By incorporating the above analysis and the ap- 
propriate screening terms in Eq. (2) we see that the 
normal-superfluid conversion is cut off at temperatures 
(or frequencies) T ~ Qa T , where Co is the bare fugacity 
of phase slips, and a~ x ~ (£/c) _1 is the UV cutoff in (2). 

Thus a finite continuous wire is effectively described 
by a chain of Josephson junctions, shunted by the 
global resistance in the chain. Phase slips now exhibit 
an interaction due to the plasma waves in the chain, 
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= , and also due to the dissipation through 
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the normal resistance, which is couple through the leads: 
| log pj . Naively, a transition will now occur when: 
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(see Fig. 4a). Less naive considerations show that the 
KTB transition, tuned by if is a cross-over for lengths 
L < c/hT, and an even stronger effect may appear due 
to the bare fugacity of phase slips being exponentially 
suppressed with K. The total resistance, however, still 
drives a Schmid transition when -Rtotai = Rq- 

Our conclusions could be easily related to the nanowire 
experiments [14-16, 32]. In long wires, we expect a SC- 
Normal crossover tuned by stiffness (as in [11, 23]), but 
in short wires, we expect a Schmid transition tuned by 
the total normal-part resistance. In Fig. 4b we recast 
the diagram of Fig. 4a for the MoGe nanowire experi- 
ments, plotting L/R tota x cx A vs. L, with L the length 
of the wire, and A its cross section area. The diago- 
nal line marks i?totai = Rq- Above it we expect T = 
superconductivity. The horizontal line marks the SC- 
normal cross-over in longer wires. This line most prob- 
ably arises from the exponential dependence of the bare 
quantum-phase-slip fugacity on thickness [18], but may 
also be associated with a KTB transition at K ~ 4, or a 
fermionic Tq suppression mechanism, which also depends 
on R^/Rq [25]. After completing the analysis described 
here, Bezryadin and coworkers measured a large number 
of short samples with L < 150nm. These show near per- 
fect fit with our prediction of a universal transition at 
-Rtotai = Rq for shorter wires [16, 32] [37]. 

The application of our simple theory to the nanowire 
experiments requires several caveats. (1) It is natural 
to associate the resistance of the nanowire devices at 
temperatures just below the SC transition of the leads, 
with the total normal-part nanowire resistance, i?totai- 
It is unclear, however, how this resistance is related to 
the normal-state resistance of the nanowires at temper- 
atures above the bulk critical temperatures for MoGe. 
(2) In addition, the origin and precise nature of nor- 
mal electrons in the wires is unknown. Possibly, phase 
fluctuations or the strong disorder stifle the proximity 
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FIG. 4: (a) Phase diagram of the two-fluid chain and wire. 
The black line describes (roughly) the SC-normal boundary 
of an infinite two-fluid chain (Fig. 2b), with r < Rq; both 
K — 2ny/ Ej I Ec and R are local quantities. When r > Rq, 
normal-superfluid conversion is cut off, the local dissipation 
becomes unimportant, and only the horizontal line applies. 
In finite continuous wires (grey line), the transition takes a 
dissipative nature when the wire is short, and occurs when 
fltotai = Rq\ longer wires have a local crossover - dashed line 
- tuned by the superconducting stiffness ~ K. (b) Phase di- 
agram for MoGe nanowires of Refs. [14-16]. The y-axis, 
L/Rtotai on A/pMoGe is proportional to the surface area, 
which is proportional to K. The x-axis is the length. The blue 
dots are insulating samples, whereas the red dots are super- 
conducting. The horizontal dashed line marks the long-wires 
cross over, while the diagonal black line marks the transition 
line Rtotai = Rq ■ 



effect and give rise to a normal part. Particularly, if 
the normal-superconductor relaxation rate is indeed sup- 
pressed, each phase slip gives rise to a long-lived pop- 
ulation of quasi-particles, as is the case in Ref. [33] 
where dissipation at phase-slip centers is investigated. 
(3) The resistance vs. temperature curves measured by 
Bezryadin on the superconducting side show sharp expo- 
nential, activated-like, decay of the resistance, contrary 
to a naive quantum phase-slip theory, where an algebraic 
dependence of the resistance on temperature is expected. 
Similarly, the wires remaining normal show a weakly in- 
sulating behavior, with a charge-gap that corresponds 
to the Coulomb-blockade of the leads [32]. These ob- 
servations do not contradict the possibility of a Schmid 
transition, and can be understood by also considering 



the effective dissipation produced by a finite density of 
phase slips, which is too large to justify the perturba- 
tive analysis pursued here. Such considerations appeared 
to describe the intermediate-coupling regime of the two- 
junction system [29]. Using the results presented here 
regarding the Schmid transition, but adding a finite den- 
sity of phase-slips, we successfully explained the sharp 
temperature dependence of short wires in Refs. [34, 35]. 

The main result of this paper is the divergence of the 
normal-superfluid relaxation time, r r = Y , in contin- 
uous uniform nanowires due to quantum fluctuations. 
Apart from the direct application of our theory to the 
MoGe nanowire experiments, this effect could be di- 
rectly investigated in meso and nanoscopic systems where 
quantum fluctuations are apparent at relatively high tem- 
peratures. Some early experiments in this direction on 
nanostructures not uniform enough, but with quantum 
fluctuations are described in Ref. [36]. In future work 
we hope to address the issues of the origin and nature 
of the normal-part in nanowires, and its interplay with 
phase-slip density, and the diverging relaxation time. 
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